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Context. Cool white dwarfs are reliable and independent stellar chronometers. The most common white dwarfs have carbon-oxygen 
, dense cores. Consequently, the cooling ages of very cool white dwarfs sensitively depend on the adopted phase diagram of the carbon- 

. oxygen binary mixture. 

Aims. A new phase diagram of dense carbon-oxygen mixtures appropriate for white dwarf interiors has been recently obtained using 
(— I ' direct molecular dynamics simulations. In this paper, we explore the consequences of this phase diagram in the evolution of cool 

white dwarfs. 

Methods. To do this we employ a detailed stellar evolutionary code and accurate initial white dwarf configurations, derived from the 
O ' full evolution of progenitor stars. We use two different phase diagrams, that of Horowitz et al. (2010), which presents an azeotrope, 

and the phase diagram of Segretain & Chabrier (1993), which is of the spindle form. 
C/3 , Results. We computed the evolution of 0.593 and 0.878 Mq white dwarf models during the crystallization phase, and we found that 

■ the energy released by carbon-oxygen phase separation is smaller when the new phase diagram of Horowitz et al. (2010) is used. 
This translates into time delays that are on average a factor ~ 2 smaller than those obtained when the phase diagram of Segretain & 
Chabrier (1993) is employed. 

, Conclusions. Our results have important implications for white dwarf cosmochronology, because the cooling ages of very old white 

I . dwarfs are different for the two phase diagrams. This may have a noticeable impact on the age determinations of very old globular 

' clusters, for which the white dwarf color-magnitude diagram provides an independent way of estimating their age. 

• Key words, stars: evolution — stars: interiors — stars: white dwarfs 
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1. Introduction is latent heat, while the second one is the release of gravitational 

^""^ . energy resulting from the changes in the carbon-oxygen profile 

■ White dwarf stars constitute the most common end-point of stel- due to the different chemical compositions of the liquid and solid 
j>! ; lar evolution — see, for instance, Althaus et al. (2010a) for a phases (Garci'a-Berro et al. 1988a, 1988b). Generally speaking, 

k> recent review — and as such ai-e valuable in constraining sev- the solid formed upon cyrstallization is richer in oxygen than the 
^ eral properties of a wide variety of stellar populations includ- liquid. As the oxygen-rich solid core grows at the center of the 
H [ ing globular and open clusters (Von Hippel & Gilmore 2000; white dwarf, the lighter carbon-rich liquid mande left behind is 
5r .Hansen et al. 2007; Winget et al. 2009; Garcia-Berro et al. efficiently redistributed by Rayleigh-Taylor instabihties (Isern et 
2010). Additionally, they can be used to place constraints on al. 1997). This process releases gravitational energy, and this ad- 
exotic elementary particles (Isern et al. 1992; Corsico et al. ditional energy source has a substantial impact in the computed 
2001 ; Isern et al. 2008) or on alternative theories of gravitation cooling times of cool white dwaifs (Segretain et al. 1994; Salai'is 
(Garci'a-Berro et al. 1995; Benvenuto et al. 2004; Garcfa-Berro et al. 1997; Montgomery et al., 1999; Salaris et al. 2000; Isern et 
et al. 2011). This is possible because we have a relatively pre- al. 2000; Renedo et al. 2010). 
cise knowledge of the main physical processes responsible for 

their evolution, although some uncertainties still persist for key Recently, Winget et al. (2009) have used theoretical fits to the 

aspects of their constitutive physics. One of the processes that observed white dwarf luminosity function of the globular cluster 

is still subject to some uncertainties is crystallization. As early NGC 6397 to provide evidence for the occurrence of crystalliza- 

recognized (Van Horn 1968), cool white dwarfs are expected to tion in deep interiors of white dwarfs, and to place constraints 

crystallize as a result of strong Coulomb interactions in their on the crystallization temperature of the carbon-oxygen mixture, 

very dense interior. Crystallization results in two additional en- Thus, it is foreseable that in the near future deep photometry of 

ergy sources that slow down the cooling process. The first source the cooling sequence of degenerate stars of nearby globular clus- 
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ters will allow us to check the accuracy of the theoretical cool- 
ing sequences, and hence will allow us to study the constitutive 
physics of matter at the high densities of cool white dwarfs. 

Since the pionering works of Stevenson (1980) and 
Mochkovitch (1983), large theoretical efforts have been paid to 
study the phase diagram of carbon-oxygen mixtures. In these 
early efforts the adopted carbon-oxygen phase diagram had a 
deep eutectic. This resulted in a high enhancement of the oxy- 
gen abundance of the solid phase. Consequently, the computed 
delays in the cooling ages were rather large. The calculations 
of Stevenson (1980) were subsequently improved by Ichimaru 
et al. (1988), who obtained a phase diagram of the azeotropic 
form, and Barrat et al. (1988), who found that the phase diagram 
was of the spindle form. Later, Ogata et al. (1993) used Monte 
Carlo simulations and the hypernetted-chain approximation to 
compute the phase diagram of the carbon-oxygen binary mix- 
ture, and obtained a phase diagram of the spindle form. Almost 
simultaneously, Segretain & Chabrier (1993) used a density- 
functional approach to obtain the phase diagram for arbitrary bi- 
nary mixtures as a function of Z\ IZ2, being Zi and Z2 the charge 
of the two chemical species. In the case of carbon-oxygen mix- 
tures, they obtained a phase diagram of the spindle type. Since in 
the case of a phase diagram of the spindle form the solid phase is 
less oxygen-enriched, the time delays computed using this type 
of phase diagrams turned out to be smaller than those previously 
computed (Segretain et al. 1994). Since then, the phase diagram 
of Segretain & Chabrier (1993) has remained a "de facto" stan- 
dard over the years. However, very recently the phase diagram of 
dense carbon-oxygen mixtures appropriate for white dwarf star 
interiors has been re-examined by Horowitz et al. (2010). This 
work was motivated by the unexpected finding of Winget et al. 
(2009) that the crystallization temperature of white dwarfs in the 
globular cluster NGC 6397 was close to that for pure carbon. 
Horowitz et al. (2010) used an approach completely different of 
those previously employed. Specifically, they used direct two- 
phase molecular dynamics simulations for the solid and liquid 
phases. Their results are in rather good agreement with those 
of Medin & Cumming (2010), and predict crystallization tem- 
peratures considerably lower than those obtained by Segretain 
& Chabrier (1993). In particular, Horowitz et al. (2010) found 
the crystallization temperature of carbon-oxygen mixtures with 
equal mass fractions to be close to that of pure carbon, thus of- 
fering a possible explanation for the puzzling result of Winget 
et al. (2009). Additionally, Horowitz et al. (2010) found that the 
shape of the carbon-oxygen phase diagram is of the azeotropic 
form, and not of the spindle type, as previously thought. This 
may have a large effect on the evolutionary ages of cool white 
dwarfs. 

In this paper we explore the implications for the evolutionary 
properties of white dwarfs of the new phase diagram of carbon- 
oxygen mixtures computed by Horowitz et al. (2010). To this 
end, we employ a detailed stellar evolutionary code and initial 
accurate white dwarf structures derived from the full evolution 
of progenitor stars. The paper is organized as follows. In Sect.|2] 
we briefly describe the phase diagram of Horowitz et al. (2010). 
We also comment on the main characteristics of our evolution- 
ary code. Particularly we also elaborate on the treatment of the 
energy sources resulting from crystallization, and the evolution- 
ary sequences used in our study. Sect. [3] is devoted to explore 
the consequences of the new phase diagram for the evolutionary 
times of white dwarfs. Finally, in Sect. 4 we summarize the main 
results of our calculations and we draw our conclusions. 
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Fig. 1. Carbon-oxygen phase diagrams used in our evolution- 
ary calculations. The crystallization temperature of the carbon- 
oxygen binary mixture in terms of the crystallization tempera- 
ture of pure carbon is shown as a function of the oxygen abun- 
dance by mass. Dashed red lines and solid blue lines correspond, 
respectively, to the phase diagrams of Segretain & Chabrier 
(1993) and Horowitz et al. (2010). For each diagram, the upper 
curve gives the crystallization temperature for a given oxygen 
abundance of the liquid, while the lower curve provides the equi- 
librium oxygen abundance of the solid at this temperature. The 
dotted line coiTesponds to the case in which no phase separation 
occurs, and the mixture is treated as the average of the chem- 
ical species. Finally, the circles denote the crystallization tem- 
perature of the carbon-oxygen mixture for the case of an oxygen 
mass abundance of 0.7, a typical value found in the evolutionary 
calculations of white dwarf progenitors. 



2. Numerical tools 

2A. Carbon-oxygen phase diagrams 

Horowitz et al. (2010) have recently determined the phase di- 
agram of dense carbon-oxygen mixtures appropriate for white 
dwarf star interiors. They argue that their results, based on di- 
rect two-phase molecular dynamics simulations, are less affected 
by small errors in the free energy difference between the liquid 
and solid phases than previous studies. Moreover, they predict 
that the crystallization temperature of the carbon-oxygen binary 
mixture is considerably lower than that resulting from the phase 
diagram of Segretain & Chabrier (1993). This can be seen Fig.[T] 
in which several phase diagrams for the carbon-oxygen mixture 
are displayed. In particular, in this figure we show both the crys- 
tallization temperature (upper curves) as a function of the oxy- 
gen abundance by mass of the liquid phase and the equilibrium 
abundances in the solid phase (lower curves) at this temperature. 
The crystallization temperature T of the carbon-oxygen mixture 
is expressed in terms of the crystallization temperature of pure 
carbon (Tc). Crystallization for pure carbon occurs when the 
Coulomb coupling parameter (F) reaches the value Fci-ys = 178.4 
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(Horowitz et al. 2010), where the Coulomb coupling parameter 
is defined as: 

being the interelectronic distance, {Z^^^) an average (by num- 
ber) over the ion charges, and Boltzmann's constant. The rest 
of the symbols have their usual meaning. In particular, the crys- 
tallization temperature for pure carbon composition results 

rc ^ 6^^^ 2.275 X 10^ ^'1 (2) 

A crys 

The shape of the phase diagram for a carbon-oxygen mix- 
ture as calculated by Horowitz et al. (2010) is displayed in Fig.[T] 
using solid blue lines. In addition, we also show (dashed red 
lines) the carbon-oxygen phase diagram of Segretain & Chabrier 
(1993). Finally, the crystallization temperature when phase sepa- 
ration is not taken into account, and the mixture is treated as the 
average of the chemical species, is represented as well (black 
dotted line). It is worth highlighting some important features of 
the phase diagrams illustrated in Fig.[r| To begin with, we note 
that the phase diagram of Horowitz et al. (2010) predicts a crys- 
tallization temperature of the carbon-oxygen mixture substan- 
tially lower than that predicted by the phase diagram of Segretain 
& Chabrier (1993), which in turn is also much lower than that 
obtained in the case in which phase separation is disregarded. 
This implies that for a white dwarf of a given mass and core com- 
position, crystallization will set in at smaller stellar luminosities 
when the phase diagram of Horowitz et al. (2010) is adopted. 
The second relevant point is that the width of the phase diagram 
of Segretain & Chabrier (1993) is considerably larger than that 
of Horowitz et al. (2010). Hence, the oxygen enhancement in the 
solid phase (with respect to the composition of the fluid phase 
from which it formed) will be substantially smaller when the 
phase diagram of Horowitz et al. (2010) is adopted. These two 
dififerences have opposite effects. On the one hand, the gravita- 
tional energy released by carbon-oxygen phase separation will 
be smaller in the case in which the phase diagram of Horowitz 
et al. (2010) is adopted. On the other, this energy is released 
at smaller luminosities. The impact of these two effects on the 
delays introduced by phase separation upon crystallization can 
only be reliably assessed using a full evolutionary code. Finally, 
there is as well another effect which is also worth noting. The 
phase diagram of Horowitz et al. (2010) presents an azeotrope 
at Xo ~ 0.3, while that of Segretain & Chabrier (1993) is ap- 
proximately of the spindle form. This means that when the fluid 
phase reaches the azeotropic composition, the solid phase has 
the same composition of the liquid. Consequently, phase separa- 
tion no longer occurs, and the subsequent evolution during the 
crystallization phase is only driven by the release of latent heat. 
Nevertheless, if the abundance of oxygen in the outer layers is 
smaller than the azeotropic one, the solid that forms is oxygen 
poor as compared with the liquid and raises until it melts, leaving 
behind a liquid that gradually aproaches to the azeotropic com- 
position. When this happens a solid with the azeotropic compo- 
sition forms. Therefore, in this case, the process of separation 
continues until the entire white dwarf core has frozen. 

2.2. Evolutionary code 

All the calculations reported here have been done using the 
LPCODE stellar evolutionary code. This code has been used to 



study different problems related to the formation and evolution 
of white dwarfs — see Althaus et al. (2010b), Renedo et al. 
(2010), and references therein for details. The only difference 
with respect to previous evolutionary calculations of cooling 
white dwarfs is the treatment of crystallization. For the results 
presented here, crystallization sets in according to the phase di- 
agram considered. That is, when no phase separation is assumed 
the crystallization temperature of the carbon-oxygen core is ob- 
tained from Eq. ([TJ, and imposing F - 180, while in all other 
cases this temperature is obtained from the corresponding phase 
diagram. 

The energy sources associated to the crystallization of the 
white dwarf core comprise the release of latent heat and grav- 
itational energy associated with changes in the carbon-oxygen 
profile induced by crystallization. In LPCODE, the inclusion of 
these two energy sources is done self-consistently and locally 
coupled to the full set of equations of stellar evolution. That is, 
the structure and evolution of white dwarfs is computed with the 
changing composition profile and with the luminosity equation 
appropriately modified to account for both the local contribution 
of energy released from the core chemical redistribution and la- 
tent heat. At each evolutionary timestep, the crystallization tem- 
perature and the change of the chemical profile resulting from 
phase separation are computed using the appropriate phase dia- 
gram. In particular, the carbon-enhanced convectively-unstable 
liquid layers overlying the crystallizing core are assumed to be 
instantaneously mixed, a reasonable assumption considering the 
long evolutionary timescales of white dwarfs (Isern et al. 1997). 
After computing the chemical compositions of both the solid and 
liquid phases the net energy released in the process is assessed as 
in Isern et al. (1997; 2000). The latent heat contribution is taken 
to be Q.llk^T per ion (Salaris et al. 2000). Both energy contri- 
butions are distributed over a small mass range around the crys- 
tallization front. We mention that the magnitude of both energy 
sources are calculated at each iteration during the convergence of 
the model — see Althaus et al. (2010c) for the numerical details. 

The input physics of the code include the equation of state 
of Segretain et al. (1994) for the high-density regime — which 
accounts for all the important contributions for both the liq- 
uid and solid phases (Althaus et al. 2007) — complemented 
with an updated version of the equation of state of Magni & 
Mazzitelli (1979) for the low-density regime. Radiative opacities 
are those of OPAL (Iglesias & Rogers 1996), including carbon- 
and oxygen-rich compositions, complemented with the low- 
temperature opacities of Alexander & Ferguson (1994), whilst 
conductive opacities are taken from Cassisi et al. (2007). During 
the white dwarf regime, the metal mass fraction Z in the en- 
velope is not assumed to be fixed. Instead, it is specified con- 
sistently according to the prediction of element diffusion. To 
account for this, we considered radiative opacities tables from 
OPAL for arbitrary metallicities. For effective temperatures less 
than 10,000 K, outer boundary conditions for the evolving mod- 
els are given by detailed non-gray model atmospheres that incor- 
porate non-ideal effects in the gas equation of state and chemi- 
cal equilibrium (based on the occupation probability formalism), 
radiative and convective transport (mixing length theory) of en- 
ergy, collision-induced absorption from molecules, and the Lya 
quasi-molecular opacity (Rohrmann et al. 2011). This provides 
detailed and accurate outer boundary conditions which are re- 
quired for a proper treatment of the evolutionary behavior of cool 
white dwarfs. 
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Fig. 2. Inner oxygen distribution (abundance by mass) for the 
0.593 Mq white dwarf model at different evolutionary stages. 
The (red) dot-dashed line corresponds to the oxygen distribution 
before the occurrence of crystallization. The final oxygen profile 
when most of the white dwarf has crystallized is also shown for 
the phase diagram of Horowitz et al. (2010) — (blue) dashed 
line — and for the phase diagram of Segretain & Chabier (1993) 
— (black) solid Une. 



2.3. Evolutionary sequences 

Instead of exploring the evolution of white dwarfs with differ- 
ent arbitrary chemical profiles for the core, we focus on evolu- 
tionary sequences that are provided by detailed calculations of 
the evolutionary history of progenitor stars. Our aim is to ex- 
plore the impact of the new phase diagram of Horowitz et al. 
(2010) on existing grids of white dwarf evolutionary sequences 
that incorporate a realistic chemical profile in the stellar interior. 
Specifically, the white dwarf initial configurations considered in 
this investigation are those obtained from the full evolution of 
progenitor stars we computed in previous studies (Renedo et 
al. 2010). In those studies progenitor stars were evolved from 
the zero age main sequence, through the thermally-pulsing and 
mass-loss phases on the asymptotic giant branch (AGB), to the 
white dwarf cooling phase. Extra-mixing episodes beyond each 
formal convective boundary were taken into account, particu- 
larly during the core helium burning stage, but not during the 
evolutionary stages corresponding to the thermally-pulsing AGB 
phase. Moreover, the outer chemical profiles of our white dwarf 
sequences are the result of element diffusion processes that lead 
to the formation of pure hydrogen envelopes — see Althaus et 
al. (2010b) for details. 

For the sake of conciseness, we concentrate on two white 
dwarf sequences, with masses 0.593 and 0.878 Mq, which are 
the result of the evolution of 1.75 and 5.0 M0 progenitors with 
metallicity Z - 0.01. The total mass of hydrogen in their en- 
velopes is 1.1 X 10"^ and 1.17 x 10"^ M0, respectively. For 
each stellar mass, we have computed the white dwarf cooling 
phase down to very low luminosities, when most of the white 
dwarf has already crystallized. Each sequence has been com- 
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Fig. 3. Same as Fig. |2]but for 0.878 Mq white dwarf model se- 
quence. 



puted considering the phase diagrams of Segretain & Chabrier 
(1993) and Horowitz et al. (2010). In the interests of compari- 
son, we have computed additional evolutionary sequences con- 
sidering the chemical profiles obtained using these two phase di- 
agrams but adopting a crystallization temperature resulting from 
imposing V - 180. In this way the effects of chemical differenti- 
ation upon crystallization can be disentangled from those result- 
ing from a smaller crystallization temperature. 



3. Evolutionary results 

We start by examining the impact of the shape of the phase dia- 
gram on the oxygen abundance distribution expected in the inte- 
rior of a crystallized white dwarf. This is illustrated in Figs.|2]and 
[3] which display the inner oxygen abundance profile at differ- 
ent evolutionary stages for the 0.593 and 0.878 Mq white dwarf 
models, respectively. In each figure, the dot-dashed line shows 
the oxygen distribution before the onset of crystallization, and 
the dashed and solid lines show, respectively, the oxygen distri- 
bution after crystallization is almost complete, when the phase 
diagrams of Horowitz et al. (2010) and Segretain & Chabrier 
(2003) are employed. As anticipated in Sect. |2] because of the 
very different shapes of both phase diagrams, we expect a dis- 
tinct oxygen distribution in the white dwarf interior by the end 
of the crystallization process. This is apparent from Figs.|2]and[3] 
where it can be seen that the final oxygen distributions predicted 
by the two phase diagrams are clearly different. Note that, for 
both masses, the phase diagram of Horowitz et al. (2010) yields 
smaller oxygen abundances in the central regions, compared 
with those obtained using the phase diagram of Segretain & 
Chabrier (1993). In particular, for the 0.593 Mq model, the crys- 
tallization process increases the oxygen abundance at the center 
by ~ 10% when the phase diagram of Horowitz et al. (2010) 
is used, and by ~ 19% when the phase diagram of Segretain & 
Chabrier (1993) is adopted. These figures turn out to be ~ 11% 
and ~ 27%, respectively, for the 0.878 Mq model white dwarf. 
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Fig. 4. Surface luminosity versus age for the 0.593 Mq white 
dwarf sequences undergoing carbon-oxygen phase separation. 
Solid and dashed curves correpond to the predictions of the 
phase diagrams of Segretain & Chabrier (1993) and Horowitz 
et al. (2010), respectively. The cooling curve in the case that 
carbon-oxygen phase separation is not considered is also shown 
as a dotted line. 



Clearly, the amount of matter redistributed by phase sepa- 
ration during crystallization for a given stellar mass is markedly 
smaller when the phase diagram of Horowitz et al. (2010) is con- 
sidered. This results in a smaller energy release from carbon- 
oxygen differentiation. Since for the 0.59 Mq model the ini- 
tial oxygen abundance is Xq ~ 0.7 when crystallization sets 
in (see Fig. |2]l, the crystallization temperature is not too dif- 
ferent for both phase diagrams (see Fig. [1]). Specifically, for 
the phase diagram of Horowitz et al. (2010) the crystallization 
temperature is ~ 1.19rc, whereas for that of Segretain et al. 
(1993) is ~ 1.35rc. Consequently, this energy is released at 
similar luminosities — log(L/LQ) ^ -3.84 for the phase di- 
agram of Horowitz et al. (2010) and log{L/LQ) ^ -3.70 for 
that of Segretain & Chabrier (1993) — and the impact on the 
white dwarf cooling times is smaller in the case in which the 
phase diagram of Horowitz et al. (2010) is adopted. The same 
occurs for the more massive white dwarf (see Fig. O. All this 
becomes clear by examining Figs.|4]and|5]which depict, respec- 
tively, the relationship between the surface luminosity and age 
for the 0.593 and 0.878 Mq white dwarf sequences that undergo 
carbon-oxygen phase separation. In each figure, solid and dashed 
curves correpond to the predictions of the phase diagrams of 
Segretain & Chabrier (1993) and Horowitz et al. (2010), respec- 
tively. In addition, the cooling curve obtained when the carbon- 
oxygen phase separation upon crystallization is neglected is 
shown with a dotted line — for this sequence, crystallization and 
the release of latent heat are assumed to occur at F = 180. Note 
that for both stellar masses, the phase diagram of Horowitz et al. 
(2010) results in white dwarf cooling times that are smaller than 
those predicted by the phase diagram of Segretain & Chabrier 
(1993). This is, as mentioned, a consequence of the smaller com- 
position changes predicted by the Horowitz et al. (2010) phase 
diagram. 



The dependence of the age delays induced by chemical dif- 
ferentiation on stellar mass for both phase diagrams is more dif- 
ficult to assess, since they depend not only on gravity but also on 
the temperature at which crystallization occurs. Since gravity is 
larger for the more massive white dwarf, the energy released by 
chemical redistribution is larger as well. However, for the more 
massive white dwarf crystallization takes place at higher stel- 
lar luminosities, and the delay in the cooling times introduced 
by chemical differentiation is smaller in this case — see, for in- 
stance, Salaris et al. (1997). In addition, the magnitude of the age 
delays are influenced by the initial chemical profile of the white 
dwarf, which is different for each stellar mass. This is more im- 
portant in the case of the phase diagram of Segretain & Chabrier 
(1993), for which the age delays strongly depend on the initial 
composition (Salaris et al. 2000), but less relevant in the case 
of the phase diagram of Horowitz et al. (2010), for which the 
predicted composition changes are smaller, see Fig.[T] 

The situation can be clarified with the help of Fig.|6l which 
illustrates the age delays resulting from carbon-oxygen phase 
separation (in percentage with respect to the case in which phase 
separation is disregarded) for the 0.593 and 0.878 Mq white 
dwarf cooling sequences (top and bottom panel, respectively), 
for both the phase diagram of Horowitz et al. (2010) — dashed 
lines — and that of Segretain & Chabrier (1993) — solid lines. 
Note that when the phase diagram of Horowitz et al. (2010) is 
employed, the age delays are negative at moderately high lumi- 
nosities, between -3.8 S log(L/L0) ^ -4.1 for the 0.593 
model. This is simply because the phase diagram of Horowitz 
et al. (2010) predicts a lower crystallization temperature than 
that obtained using Eq. ([1]) with F = 180. Hence, crystalliza- 
tion occurs at lower stellar luminosities when the phase dia- 
gram of Horowitz et al. (2010) is employed, with the conse- 
quence that the coiTesponding cooling sequence has initially 
shorter cooling ages. It is also worth emphasizing that when 
the phase diagram of Horowitz et al. (2010) is adopted, the 
bulk of the delay in the cooling ages of an otherwise typical 
white dwarf of 0.593 Mq occurs at luminosities ranging from 

log(L/LQ) 4.2 to log(L/LQ) 4.5, whereas for the phase 

diagram of Segretain & Chabrier (1993) this occurs at luminosi- 
ties between log(L/L0) 4.0 and -4.4, when a sizable part of 

the core is already crystallized. But the most relevant aspect that 
Fig. |6]puts forward is that when the phase diagram of Horowitz 
et al. (2010) is adopted, carbon-oxygen phase separation is less 
relevant for the evolution of white dwarfs. In particular, for the 
0.593 Mq at log(L/L0) ^ -4.5, the age delay is only ~ 0.5 Gyr 
when the phase diagram of Horowitz et al. (2010) is adopted, 
whilst in the case in which that of Segretain & Chabrier (1993) 
is used the delay amounts to ~ 1.3 Gyr. 

We also note that for the 0.593 Mq white dwarf cooling se- 
quence, the phase diagram of Horowitz et al. (2010) predicts 
ages that are up to ^ 8% shorter at log(L/Lo) = ~4.0 when 
compared with the ages derived using the phase diagram of 
Segretain & Chabrier (1993). The differences are larger for the 
0.878 Mq white dwarf sequence, reaching up to ^ 17% also at 
log(L/L0) = -4.0. In this case the differences are larger because 
the smaller initial oxygen abundance translates into a larger rel- 
ative enrichment of the solid phase when the phase diagram of 
Segretain & Chabrier (1993) is used — but not in the Horowitz et 
al. (2010) phase diagram. This, in turn, results in a larger energy 
release, and consequently in larger delays. This merely reflects 
the fact that age delays resulting from carbon-oxygen separation 
are less sensitive to the initial chemical profile when using the 
Horowitz et al. (2010) phase diagram. 
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Age (Gyr) 

Fig. 5. Same as Fig. |4]but for the 0.878 white dwarf se- 
quences. 




Fig. 6. Age delays, in percentage with respect to the case in 
which phase separation is not considered, due to carbon-oxygen 
phase separation according to the phase diagram of Horowitz et 
al. (2010) and Segretain & Chabrier (1993), dashed and soUd 
lines, respectively. Upper and bottom panel correspond to the 
0.593 and 0.878 Mq white dwarf sequences, respectively. 



To isolate the effect on the cooling ages of the shape of the 
phase diagram from the different crystallization temperature of 
both phase diagrams, in Fig. Q we compare the cooling curves 
for several 0.593 Mq white dwarf sequences. In addition to the 
cooling curves considered in Fig. 31 we also plot the cooling 
curves computed for both phase diagrams when the crystalliza- 
tion temperature is given by the relation F - 180, rather than 



phase diagrams and white dwarf evolution 




Age (Gyr) 

Fig. 7. Same as Fig. |4]but including the cooling curves (thin 
lines) that result from imposing that crystalUzation occurs at 
r= 180. 



being obtained using the upper curve of the phase diagram. 
These cooling curves are shown using thin lines. The results 
deserve some comments. As already noted, the phase diagram 
of Horowitz et al. (2010) predicts lower crystallization tem- 
peratures (and luminosities) than that of Segretain & Chabrier 
(1993), and even lower than the crystallization temperature when 
no phase separation is considered. Specifically, for the 0.593 Mq 
sequence, crystallization starts at log{L/LQ) = -3.70 when the 
phase diagram of Segretain & Chabrier (1993) is used, and at 
log(L/L0) = -3.84 when that of Horowitz et al. (2010) is 
adopted, while when no phase separation is considered, crystal- 
lization occurs at log(L/L0) = -3.68. This explains the initially 
larger cooling ages in the latter sequence. But at smaller lumi- 
nosities, the trend is reversed. We thus conclude that it is the 
shape of the phase diagram that is the most relevant factor in- 
fluencing the delays in the ages of cool white dwarfs, and not 
the specific value of the crystallization temperature. Finally, it is 
worth mentioning as well that the value of F at the onset of crys- 
tallization is 186 and 210 for the phase diagrams of Segretain & 
Chabrier (1993) and Horowitz et al. (2010). These values turn 
out to be somewhat higher, 190 and 218, respectively, for the 
0.878 Mq sequence, because of the initially lower oxygen abun- 
dance of this sequence. 

4. Conclusions 

Winget et al. (2009) suggested that the crystallization tempera- 
ture of carbon-oxygen white dwarf cores in the globular cluster 
NGC 6397 is close to that of pure carbon. This unexpected re- 
sult prompted Horowitz et al. (2010) to determine a new phase 
diagram for carbon-oxygen mixtures using direct molecular dy- 
namics simulations for the solid and liquid phases. They found 
crystallization temperatures considerably lower than those given 
by the most usually adopted prescription, which is obtained im- 
posing F = 180, and neglecting carbon-oxygen phase separation. 
In particular, Horowitz et al. (2010) found that the crystallization 
temperature for carbon-oxygen mixtures with equal mass frac- 
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tions of carbon and oxygen to be close to that of pure carbon, 
thus offering an explanation for the puzzling result of Winget et 
al. (2009). They also found that the shape of the phase diagram 
of binary carbon-oxygen mixtures is of the azeotropic form, in- 
stead of the spindle phase diagram (Segretain & Chabrier 1993) 
previously employed in the most accurate calculations of cooling 
white dwarfs (Renedo et al. 2010) available so far. The core fea- 
ture of this paper has been precisely to explore the consequences 
for white dwarf evolution of this new phase diagram. To this end, 
we used the LPCODE stellar evolutionary code, and we computed 
several cooling sequences for white dwarfs of masses 0.593 and 
0.878 Mq. The initial white dwarf configurations were extracted 
from the full and detailed evolution of progenitor stars we com- 
puted in previous studies, which also provided realistic initial 
chemical profiles. 

The lower crystallization temperature predicted by the 
Horowitz et al. (2010) phase diagram means that the onset 
of crystallization in white dwarfs occurs at stellar luminosities 
smaller than those predicted by the phase diagram of Segretain 
& Chabrier (1993). For a typical white dwarf of mass ~ 0.6 Mq, 
we find that the phase diagram of Horowitz et al. (2010) predicts 
crystallization to occur at log(L/LQ) = -3.84, while the lumi- 
nosity at which crystalUzation sets in when the phase diagram of 
Segretain & Chabrier (1993) is used is log(L/L0) = -3.70, and 
when no phase separation is considered the white dwarf core 
crystallizes at log(L/LQ) = -3.68. Additionally, for this new 
phase diagram the value of Coulomb coupling parameter at the 
onset of crystallization for a 0.6 Mq white dwarf is T ^ 210. 

The amount of matter that is redistributed by phase separa- 
tion during crystallization is notably smaller in the new phase 
diagram than in previous calculations of this kind. Hence, we 
find that carbon-oxygen phase separation becomes less relevant 
for white dwarf evolution when this phase diagram is adopted. 
At the luminosities for which a large fraction of the white dwarf 
mass has crystallized, we find age delays due to carbon-oxygen 
phase separation that are on average a factor ~ 2.5 smaller than 
the delays obtained using the phase diagram of Segretain & 
Chabrier (1993). Another interesting feature of the new phase di- 
agram of Horowitz et al. (2010) is that composition changes are 
less sensitive to the initial chemical profile of the white dwarf. 
This is a relevant point, since it means that the magnitude of the 
age delays induced by carbon-oxygen phase separation will be 
less aff'ected by current uncertainties in the initial carbon-oxygen 
composition of white dwarfs. 

Our results have implications for the age determinations 
of stellar populations using the white dwarf cooling sequence, 
which should be investigated in subsequent works. This is par- 
ticularly relevant for the well-studied, old, metal-poor globular 
cluster NGC 6397, which has been imaged down to very faint 
luminosities and for which a reliable color-magnitude diagram 
and a white dwarf luminosity function have been derived. As 
mentioned, this cluster has been recently used to constrain the 
properties of crystallization in the deep interior of cool white 
dwarfs (Winget et al. 2009). In particular, Winget et al. (2010) 
have shown that the observed white dwarf luminosity function in 
NGC 6397 seems to be consistent with pure carbon core white 
dwarfs crystallizing at F « 178, or, alternatively, carbon-oxygen 
core white dwarfs crystallizing at F values larger than 178, the 
theoretical value for a one component plasma. This finding is 
in line with the predicions of molecular dynamics simulations 
of Horowitz et al. (2010). These results together with the ones 
reported in this paper call for the need of studying the cooling 
sequence of crystallizing white dwarfs in other old stellar clus- 



ters on the basis of Horowitz et al. (2010) phase diagram. Work 
in this direction is in progress. 
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